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VLF/LF  Reflectivity  of  the  Polar  Ionosphere 
2 January— 30  April  1977 


1.  INTRODUCTION 

This  paper  provides  a summary  of  high  latitude  ionospheric  reflectivity,  as 
observed  by  the  USAF's  high  resolution  VLF/LF  ionosounder  operating  in  northern 
Greenland.  **  2 As  shown  in  Figure  1,  the  transmitter  is  located  at  Thule  Air  Base, 
Greenland  (76°  33’N.  Lat. , 68°  40'W.  Long.),  and  the  receiving  site  is  106  km 
north  at  the  Danish  Meteorological  Institute's  Ionospheric  Observatory  in  Qanaq, 
Greenland  (77°  24'N.  Lat.,  69°  20'W.  Long..  Geomagnetic  Lat.  89°  06'N. ).  The 
ionosounding  transmissions  consist  of  a series  of  extremely  short  (approximately 
100  psec)  VLF  pulses,  precisely  controlled  in  time,  and  radiated  from  a 130  meter 
vertical  antenna.  At  the  receiving  site,  orthogonal  loop  antennas  are  used  to  sepa- 
rate the  two  polarization  components  of  the  ionospherically  reflected  skywave  signal. 
One  antenna,  oriented  in  the  plane  of  propagation,  is  used  to  sense  the  groundwave 
and  the  "parallel"  component  of  the  downcoming  skywave.  The  second  loop,  nulled 
on  the  groundwave,  senses  the  "perpendicular"  skywave  component.  The  signal  from 
each  of  the  antennas  is  digitally  averaged  to  improve  the  signal -to -noise  ratio  of  the 
individual  received  waveforms  before  they  are  recorded  on  magnetic  tape.  An 
(Received  for  publication  28  July  1977) 

1.  Lewis.  E.A.,  Rasmussen,  J.E.,  and  Kossey,  P.A.  (1973)  Measurements  of 
ionospheric  reflectivity  from  6 to  35  kHz,  J,  Geophys.  Res.  7Jhl9. 


2.  Kossey,  P.A.,  Rasmussen,  J.  E. , and  Lewis,  E.  A.  ( 1974)  VLF  pulse  iono- 
sounder measurements  of  the  reflection  properties  of  the  lower  ionosphere, 
Akademie  Verlaq,  COSPAR,  July. 


example  of  the  observed  waveforms  is  given  in  Figure  2,  where  the  parallel  wave- 
form (Figure  2a)  consists  of  (1)  a groundwave  propagated  pulse,  (2)  a quiet  interval 
containing  low  level,  off  path  groundwave  reflections,  followed  by  (3)  the  first-hop 
skywave  signal.  The  perpendicular  waveform  is  shown  in  Figure  2b. 

Ionospheric  reflection  parameters  are  derived  by  computer  (AFGL's  CDC 
6600)  processing  of  the  ground  and  ionospherically  reflected  waveforms  with 
allowance  made  for  factors  such  as  ground  conductivity  and  antenna  patterns  (see 
Section  4). 

Although  the  data  are  recorded  about  once  per  minute,  for  this  paper  the  wave- 
forms are  averaged  into  two-hour  time  blocks  with  the  exception  of  the  three-dimen- 
sional waveform  presentations  (Section  2.2).  The  resulting  information  is  presented 


in  a weekly  format  (Figures  3 through 
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Figure  2.  Example  of  the  Observed  Waveforms 

2.  OBSERVED  WAVEFORMS 

2.1  Weekly  Example  of  Individual  Waveforms 

In  part  A of  Figures  3 fhrough  19  a set  of  averaged  parallel  and  perpendicular 
waveforms  is  presented  for  the  time  block  centered  near  local  noon  of  the  indicated 
day.  Each  of  these  waveforms  is  comprised  of  256  digitally  averai  rd  p'unts  spaced 
2 p sec  apart.  In  part  B of  the  figures,  the  groundwave  Fourier  amplitudes  are  shown 
as  a function  of  frequency.  Although  the  data  presented  in  parts  C through  E of  the 
figures  are  generally  limited  to  frequencies  in  the  first,  or  principal,  lobe  of  the 
spectrum,  information  at  higher  frequencies  can  be  used  when  sufficient  signal -to- 
noise  conditions  exist.  There  is,  however,  a frequency  range  around  each  spectral 
null  where  insufficient  signal  exists  for  measurements. 

2.2  Three-Dimensional  Waveform  I’resentation 

A three-dimensional  display  of  the  recorded  p waveforms  covering  each  weekly 
period  is  shown  in  Part  R of  each  figure  and  the  corresponding  ]_waveforms  are 
shown  in  Part  S.  For  these  plots  the  data  has  been  averaged  into  fifteen  minute 
time  blocks. 
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3.  REFLECTION  HEIGHTS 

The  group  mirror  height  (GMH)  of  reflection  was  obtained  by  determining  the 
group  delay  of  the  skywave  relative  to  the  groundwave  and  attributing  this  time 
difference,  by  simple  geometry  (assuming  a sharply  bounded  mirror-like  ionos- 
phere) to  a difference  in  propagation  distance.  As  discussed  in  Lewis  et  a I,  * the 
group  delay  can  be  defined  as  the  rate  of  change  of  phase  with  frequency.  F'or  the 
GMH  data  presented  in  this  paper,  a finite  frequency  difference  of  1.0  kHz  was 
ised,  and  the  corresponding  phase  difference  as  a function  of  frequency  for  the 
groundwave  and  both  skywave  signals  was  obtained  by  Fourier  analysis  of  the 
respective  pulses.  The  GMH  calculations  took  into  account  ground  conductivity 
(10  mho/meter  is  assumed),  and  the  corrections  of  Wail  and  Howe'  were  applied. 

Group  mirror  heights  are  plotted  as  a function  of  frequency  in  parts  C and  D 
of  Figures  3 through  19,  as  obtained  from  the  parallel  and  perpendicular  wave- 
forms, respectively.  The  GMH's  are  also  presented  as  a function  of  time-of-day 
for  the  average  frequency  of  16.5  kHz  in  figure  parts  E and  I.  The  parallel  GMH's 
in  part  E are  shown  along  with  an  average  reflection  height  for  reference  purposes. 
Each  point  of  the  reference  height  is  a weekly  average,  by  time  block,  for  the 
7 -day  period  indicated.  The  corresponding  perpendicular  GMH's,  part  I of  the  fig- 
ures, are  also  shown  with  the  weekly  average  for  comparison.  Part  G gives  the 
average,  by  time  block,  for  the  daily  parallel  GMH  data  of  part  E,  and  part  K gives 
the  corresponding  perpendicular  GMH  averages  from  the  daily  data  of  part  I. 

4.  REFLECTION  COEFFICIENTS 

Assuming  that  the  ionosphere  acts  as  a "mirror"  at  the  GMH,  plane  wave  re- 

4 

flection  coefficients  were  obtained  by  comparing  the  ratio  of  the  skywave  Fourier 
amplitude  at  a specific  frequency  to  that  of  the  groundwave,  taking  into  account  the 
antenna  patterns,  wave  spreading,  earth  curvature,  ground  conductivity,  path 
lengths,  and  antenna  patterns  including  ground  image  effects. 

The  reflection  coefficient  ||R||  was  obtained  from  analysis  of  the  parallel  sky- 
wave  component  and  is  plotted  as  a function  of  frequency  in  part  C of  Figures  3 
through  19.  The  ||R||  coefficient  for  16  kHz  is  plotted  as  a function  of  time-of-day 
in  part  F along  with  the  average  of  the  indicated  week  for  reference  purposes. 

From  the  perpendicular  skywave  pulse,  the  coefficient  ||Rj^was  obtained  andappears 

3.  Wait,  J.R.,  and  Howe,  H.  H.  (1956)  Amplitude  and  Phase  Curves  for  Ground- 

Wave  Propagation  in  the  Band  200  Cycles  per  Second  to  500  Kilocycles,  Nat. 
Bur.  Stand.  U.  S.  Circ.  No.  574.  * 

4.  Budden,  K.  G.  (1961)  Radio  Waves  in  the  Ionosphere,  p.  85,  Cambridge 

University  Press,  London. 


k . I 
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as  a function  of  frequency  in  part  D.  The  16  kHz  ||  Rj_  is  shown  along  with  its 
reference  in  part  J.  Parts  H and  L present  the  average,  by  time  block,  of  the 
daily  ||  R||  and  ||  R j_  data  presented  in  parts  F and  J,  respectively. 

For  certain  coefficient  data  points,  plotted  as  asterisks  (*),  the  reflection 
coefficient  appears  without  a corresponding  GMH.  For  these  particular  data,  only 
the  skywave-groundwave  ratios  could  be  obtained  as  the  skywaves  were  too  weak 
to  provide  reliable  group  delay  information.  The  reflection  coefficients  were  there- 
fore estimated  using  a nominal  GMH  of  80  km  in  the  calculations.  These  estimated 
coefficient  values  are  included  in  the  averages  presented  in  parts  H and  L,  but  the 
assumed  heights  are  not  used  in  the  GMH  averages  shown  in  parts  G and  K. 


5.  SUPPLEMENTARY  INFORMATION 

For  purposes  of  comparison  and  interpretation,  certain  supplementary  data  are 
presented  in  Figures  3 through  19.  Part  M of  the  figures  shows  the  magnitude  of 
the  horizontal  component  of  the  polar  magnetic  field  observed  with  a three-axis 
fluxgate  magnetometer,  and  part  N presents  30-MHz  riometer  data,  an  indicator  of 
D-region  particle  precipitation.  Figure  parts  P and  Q give  the  received  VLF  phase 
and  amplitude  from  the  18.  6 kHz  station  NLK  (transmitter  location:  Tim  Creek, 
Washington),  as  observed  at  Thule  AB  over  a 3900-km  propagation  path.  These 
supplementary  data  were  recorded  at  30 -sec  intervals  by  RADC/ETEE  at  Thule  AB; 
the  curves  represent  the  average  of  10 -min  periods.  The  solar  zenith  angle  is  given 
in  part  O for  the  indicated  mid-week  date. 


6.  ADDITIONAL  COMMENTS 

It  is  noted  that  the  transmitter  power  was  occasionally  unstable  during  the 
period  between  day  93  and  day  103  as  indicated  in  parts  R and  S of  Figures  17  and 
18.  As  the  ionospheric  reflection  height  and  reflection  coefficient  calculations  are 
based  on  comparisons  between  the  groundwave  and  skywaves,  which  are  equally 
influenced  by  the  power  variation,  these  data  remain  unaffected  through  this  period. 

This  report  is  one  of  a series.  ^ ^ Comments  and  suggestions  for 

improving  its  usefulness  should  be  addressed  to  the  VLF/ULF  Techniques  Branch 
(ETEE),  Electromagnetic  Sciences  Division,  Deputy  for  Electronic  Technology 
(RADC/ETEE),  Hanscom  AFB,  Massachusetts  01731. 

(Because  of  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  Reference  Page  No.  81,  for  References  5 through  11. 
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Figure  3.  VLF/LF  Renectivity  Data  for  the  Polar  Ionosphere,  DAY  2 (2  Jar.)  - DAY  8 (8  Jan)  1977 


Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  2 (2  .Jan)  — DAY  8 (8  Jan)  1977  (Cont) 


Figure  3.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere 
DAY  8 (8  .Ian)  1977  (Cont).  Part  R,  ||Waveform  Display 


Figure  4.  VLF/'LF  Reflectivity  Data  for  the  Polar  Ionosphere.  DAY  9 (9  Jan)  - DAY  15  (15  Jan)  1977 


Figure  4.  VLF/  LF  Reflectivity  Data  for  the  Polar  Ionosphere.  DAY  1 - .15  1077  (Co nt) 


Figure  4.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere 
DAY  15  (15  Jan)  1977  (Cont).  Part  R,  ||  Waveform  Display 
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Figure  4.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  9 (9  Jan) 
DAY  15  (15  Jan)  1977  (Cont).  Part  S,  J_Waveform  Display 
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Figure  5.  VDF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  16  (16  Jan)  — DAY  22  (22  Jan)  1977 


Figure  5.  VLF/LF  Reflectivity  Data  for  the  Polar  ionosphere 
DAY  22  (22  Jan)  1977  (Cont).  Part  R,  || Waveform  Display 


Figure  6.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  23  (23  !an)  — DAY  29  (29  Jan)  1977 
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Figure  6.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  23  (23  Janl 
DA^  29  (29  ,Tan)  1977  (Cont).  Part  R,  ||Waveform  Display 


Figure  7.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  30  (30  Jan)  — DAY  36  (5  Feb)  1977 


Figure  7.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  39  (30  Jan)  — DAY'  36  (5  Feb)  1977  (Cont) 


Figure  7.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere 
DAY  36  (5  Feb)  1977  (Cont).  Part  S,  J_Waveform  Display 


Figure  8.  VLF/LF  Reflectivity  Data  for  the'Polar  Ionosphere,  DAY  37  (6  Feb)  - DAY  43  (12  Feb)  1977 


Figure  8.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  37  (6  Feb)  — DAY  43  (12  Feb)  1977  (Cont) 


Figure  8.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere 
DAY  43  (12  Feb)  1977  (Cont).  Part  R,  || Waveform  Display 
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Figure  10.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  51  (20  Feb)  — DAY  57  (26  Feb)  1977 


Figure  10.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  51  (20  Feb)  — DAY  57  (26  Feb(  1077  (Cont) 


Figure  10.  VLF/l.F  Reflectivity  Data  for  the  Polar  Ionosphere 
DAY  57  (26  Feb)  1977  (Cont).  Part  R,  ||  Waveform  Display 


Figure  10.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  51  (20  Feb) 
DAY  57  (26  Feb)  1977  (Cont).  Part  S,  J_Waveform  Display 


Figure  11.  VLF/LF  Reflectivity  Data  Tor  the  Polar  Ionosphere 
DAY  64  (5  Mar)  1977  (Cont).  Part  S,J_Waveform  Display 
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Figure  12.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  65  (6  Alar)  - DAY  7 1 (12  Mar)  1977  (Cont) 
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Reflectivity  Data  for  the  Polar  Ionosphere 
(Cont).  Part  R,  || Waveform  Display 


Figure  14.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere 
DAY  85  (26  Mar)  1977  (Cont).  Part  R,  ||  Waveform  Display 


Figure  15.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  86  (27  Mar)  — DAY  92  (2  Apr)  1977 


Figure  15.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  86  (27  Mar)  — DAY  92  (2  Apr)  1977  (Cont) 


Reflectivity  Data  for  the  Polar  Ionosphere 
(Cont).  Part  R,  ||  Waveform  Display 


Figure  16.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  93  (3  Apr)  - DAY  99  (9  Apr)  1977 
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Figure  16.  VI.F/LF  Reflectivity  Data  for  the  Polar  Ionosphere 
DAY  99  (9  Apr)  1977  (Cont).  Part  R,  J|  Waveform  Display 


^««N"MiTT£R  POWER  N 'AR.E- 


Figure  16.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere 
DAY  99  (9  Apr)  1977  (Cont).  Part  S,  J_Waveform  Display 


Figure  17.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  100  (10  Apr)  — DAY  106  (16  Apr)  1977  (Cont) 
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Figure  18.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  107  (17  Apr)  — DAY  113  (23  Apr)  1977 


Figure  18.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere 
DAY  113  (23  Apr)  1077  (Conti.  Part  R.  II  Waveform  Display 


Figure  19.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  114  (24  Apr)  - DAY  120  (30  Apr)  1977 


Figure  19.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere,  DAY  114  (24  Apr)  — DAY  120  (30  Apr)  1977  (Cont) 


Figure  10.  VLF/LF  Reflectivity  Data  for  the  Polar  Ionosphere 
DA\  120  (30  Apr)  1077  (Cont).  Part  R,  ||  Waveform  Display 
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METRIC  SYSTEM 


BASE  UNITS: 


Quantity 

Unit 

SI  Symbol 

Formula 

length 

metre 

m 

mass 

kilogram 

kg 

.. 

time 

second 

s 

electric  current 

ampere 

A 

thermodynamic  temperature 

kelvin 

K 

amount  of  substance 

mole 

mol 

luminous  intensity 

candela 

cd 

SUPPLEMENTARY  UNITS: 

plane  angle 

radian 

rad 

solid  angle 

steradian 

sr 

DERIVED  UNITS: 

Acceleration 

metre  per  second  squared 

m/s 

activity  (of  a radioactive  source) 

disintegration  per  second 

(disintegration)^ 

angular  acceleration 

radian  per  second  squared 

rad/s 

angular  velocity 

radian  per  second 

rad/s 

area 

square  metre 

m 

density 

kilogram  per  cubic  metre 

kg/m 

electric  capacitance 

farad 

F 

A-s/V 

electrical  conductance 

siemens 

S 

A/V 

electric  field  strength 

volt  per  metre 

V'm 

electric  inductance 

henry 

II 

Vs/A 

electric  potential  difference 

volt 

V 

W/A 

electric  resistance 

ohm 

VIA 

electromotive  force 

volt 

V 

W/A 

energy 

joule 

1 

N-m 

entropy 

joule  per  kelvin 

J/K 

force 

newton 

N 

kg-m/s 

frequency 

hertz 

Hz 

( cycle I/s 

illuminance 

lux 

lx 

Im/m 

luminance 

candela  per  square  metre 

cd/m 

luminous  flux 

lumen 

Im 

cd-sr 

magnetic  field  strength 

ampere  per  metre 

A/m 

magnetic  flux 

weber 

Wb 

V-s 

magnetic  flux  density 

tesla 

T 

Wb/m 

magnetomotive  force 

ampere 

A 

power 

watt 

W 

!'» 

pressure 

pascal 

Pa 

N/m 

quantity  of  electricity 

coulomb 

C 

A-s 

quantity  of  heat 

joule 

1 

N-m 

radiant  intensity 

watt  per  steradian 

W/sr 

specific  heat 

loule  per  kilogram-kelvin 

1/kg-K 

stress 

pascal 

Pa 

N/m 

thermal  conductivity 

watt  per  metre-kelvin 

W/m-K 

velocity 

metre  per  second 

m/s 

viscosity,  dynamic 

pascal -second 

Pa-s 

viscosity,  kinematic 

square  metre  per  second 

... 

m/s 

voltage 

volt 

V 

W/A 

volume 

cubic  motre 

m 

wavenumber 

reciprocal  metre 

(wave)/m 

work 

joule 

1 

N-m 

SI  PREFIXES: 

Multiplic 

ation  Factors 

Prefix 

SI  Symbol 

1 000  000  000  000  * W’1 

t era 

T 

1 000  000  000  = 10* 

gig* 

G 

1 000  000  = 10* 

mega 

M 

1 000  = 10' 

kilo 

k 

100  = 10J 

hecto* 

h 

10  = 10' 

deka* 

da 

0 1 - io-> 

decl* 

d 

001  - io-> 

centl* 

c 

0001  = io-» 

mill! 

m 

0 000  001  * 10_* 

micro 

M 

0 000  000  001  = 10-* 

nano 

n 

0 000  000  000  001  = 10- 

pico 

p 

0 000  000  000  000  001  » 10-'* 

femto 

f 

0 000  000  000  000  000  001  - 10-'§ 

atto 

a 

* To  be  avoided  where  possible 


